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Motivated by recent measurements on CuGeC>3 the spectrum of magnetic excitations of an antifer- 
romagnetic S — \ chain with alternating coupling strength is investigated. Wave vector dependent 
magnons and a continuum with square root behavior at the band edges are found. The spectral 
density of the continua is calculated. Spin rotation symmetry fixes the gap of the continuum to 
be twice the elementary magnon gap. This is in excellent agreement with experimental results. In 
addition, the existence of bound states of two magnons is predicted: below the continuum a singlet 
and a triplet, above the continuum an "anti-bound" quintuplet. The results are based on field 
theoretic arguments, RPA calculations, and consideration of the limit of strong alternation. 
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One-dimensional quantum antiferromagnets can ex- 
hibit the interesting phenomenon of the spin-Peierls 
transition, i.e. a spontaneous dimerization of the lattice 
and the exchange constant. The recent discovery 0-^] 
of a spin-Peierls transition in the inorganic compound 
CuGeOa for the first time made neutron scattering in- 
vestigations in the ordered state possible. These provide 
interesting information on the dynamic spin correlations 
and therefore on the magnetic excitations. 

The neutron scattering experiments exhibit a number 
of interesting features: the magnon line shape asymme- 
try observed experimentally and the interpretation 
of the behavior of the equal-time structure factor in- 
dicate quite clearly the existence of a continuum besides 
a dominant magnon peak. Its weight is roughly equiva- 
lent to the one in the main peak |5| . Very recent neutron 
scattering measurements with high resolution |?J show 
a "double gap" structure: at k — tt one finds beside 
a magnon peak at the gap u) — A another gap which 
separates the magnon from the band edge of the contin- 
uum. It turns out that this second gap equals within 
experimental accuracy the first one, i.e. the onset of the 
continuum is at 2A. Furthermore, there is presently an 
active debate on the importance of next-nearest neighbor 
couplings PH10I1. Evidence for a finite competing next- 
nearest neighbor coupling stems from fits of the magnetic 
susceptibility HH x(^) an( i from a fit of the equal-time 
structure factor. 

These points motivate the present work. Its aim is 
to clarify theoretically the magnetic excitation spectrum 
in the dimerized, spin-Peierls phase below the transition 
temperature Tgp ~ 14. 3K in the dimerized phase. 

In the dimerized phase it is reasonable to adopt an adi- 
abatic approach toward the coupling between the mag- 
netic and the phonon subsystem p|-|ri]]. This is possi- 
ble since the phonon subsystem is three dimensional and 
hence its fluctuations are weak. This holds in particu- 
lar well below the transition temperature. It seems fur- 
thermore legitimate to neglect to first approximation the 
interchain coupling ||,||]. 

For the above reasons we will focus on the following 
Hamiltonian 
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where the Si are spin-1/2 operators. Beside nearest- 
neighbor coupling J > 0, next-nearest neighbor coupling 
aJ has been included. The values of a used so far range 
between 0.24 and 0.36 |p|-|l0|. We now go to a fermionic 
representation of the physics of the model via the Jordan- 
Wigner transformation Il2) to fermionic degrees of free- 
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For a/0 next-nearest neighbor terms are also present. 
For 6 = the properties of this model are well- 
understood, mainly thanks to Bethe's exact solution 
]l3| p^[ . In particular, there are no spin-wave excita- 
tions, but the low-lying excitations rather form a spinon 
continuum with lower bound (ttJ/2) sinfc. 

We now note that in the dimerized phase the momen- 
tum k is not conserved: any mode at k is coupled to the 
mode k + tt. If lo belongs to the spectrum at k, it does 
so at k + 7r as well. Together with the inversion symme- 
try k <-> — k one obtains thus reflection symmetry about 
k = tt/2. This symmetry is not taken into account in pre- 
vious results [po| , pT[ . The reflection symmetry, however, 
only refers to the energies, not to the spectral weights. 

We first consider the case of weak dimerization, 5 <C 1. 
Then a continuum approximation for the fermionic model 
(^|) is appropriate, and one obtains the Hamiltonian of the 
so-called massive Thirring model, for which there exists 
an exact Bethe ansatz solution |l7]] . The model describes 
relativistic fermions with "light velocity" v = ttJ/2 and 
a bare mass mo cx 6. In addition there is a fermion- 
fermion interaction of strength g. The known scaling re- 
lation p"8| ] for the magnetic gap (the renormalized mass 
in fermionic language), A cx S 2 ^ 3 (up to logarithmic 
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corrections) uniquely fixes g = 2v. The exact solu- 
tion provides valuable information about the spectrum: 
first there are single-particle and single-hole excitations 
which in the spin language correspond to states with 
AM Z = ±1, with excitation energy \/ A 2 + v 2 k 2 . From 
the fermionic standpoint it is then clear that there is a 
particle-hole continuum (AM Z = 0) with minimum en- 
ergy 2A. However, due to the presence of interactions in 
the model, there is a well-defined exciton-like particle- 
hole bound state below this continuum jlT]], at energy 
A when g — 2v. The particle, hole, and exciton thus 
form the low-energy degenerate triplet excitation (the 
magnon) expected because of spin rotation invariance of 
the original lattice model Q . An applied magnetic field 
would split this triplet into its three components. In 
addition, there is another excitonic bound state at en- 
ergy \/3A/2 Jl7| . This, however, has no counterpart in 
the AM Z = ±1 sectors, and therefore is a singlet exci- 
tation. Contrary to what was proposed in ref. jll[] it is 
thus not possible to continuously connect (in fc-space) 
the two bound state excitations. We note that in a 
bosonized picture of the fermionic theory |T5| , |l6| ] , the par- 
ticle, hole, and exciton correspond to soliton, antisoliton 
and breather modes. Also, including a finite a affects the 
values of v and A quantitatively, but does not affect the 
general structure of the spectrum which is protected by 
spin rotation invariance. 

From the elementary triplet at energy A it is possible 
to create multiple excitation continua. Exciting two of 
the triplet modes, nine different types of states can be 
created which form states with total spin 5 = 0, 1, and 2 
and minimal energy 2 A. Only the 5 = 1 state is of course 
accessible in a low-temperature neutron scattering exper- 
iment. It is interesting to consider what happens when 
the dimerization and therefore the gap go to zero. Then 
n-magnon continua, with lower edge at nA, become close 
to the 2-magnon continuum and they all coincide when 
A — > 0, giving rise to the well-known [ [l4| spinon contin- 
uum of the homogeneous spin-1/2 antiferromagnet. One 
should however notice that the elementary constituents 
of the continua are 5 = 1 objects when 5^0, whereas 
they are 5 = 1/2 spinons (or delocalized domain walls) 
for 5 = 0. In fact, from the bosonic representation [ |15|Jl6| 
it is clear that a nonzero S acts as a linear confining po- 
tential between two spinons, i.e. spinons can not exist as 
elementary excitations when 6^=0. 

From the above discussion the experimental observa- 
tion of a 5 = 1 continuum at twice the energy of the 
elementary magnon [jf) seems very natural. To inves- 
tigate the spectral shape of the continuum further we 
undertook a RPA study at zero temperature of the lon- 
gitudinal spin correlation function (the density correla- 
tion function in fermionic language), starting from eq.(|^). 
The general picture is the same as that of the continuum 
limit: there is a gap A in the single-particle (-hole) exci- 
tation, a particle-hole continuum above 2A, and an ex- 
citonic bound state within this gap. The existence of the 
gap for S > helps the approximation by providing an 



infrared cutoff, i.e. by suppressing low-energy particle- 
hole excitations. For isolated dimers (5 — > 1, a — 0), the 
RPA is exact. Yet in general the approximation breaks 
the spin rotational invariance since the S z component is 
dealt with differently than the S x and S v components. 
In a heuristic approach, this can be remedied by choos- 
ing an effective interaction constant J — > sJ such that 
the exciton (magnon) is in the middle of the two-particle 
continuum gap. Since the static Fock diagram renormal- 
izes the hopping one works in addition with an effective 
energy scale J e ff and an effective alternation 5 e ff. They 
are both enhanced in comparison to the bare ones. 
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FIG. 1. Left: Gap A vs. alternation S compared to pertur- 
bative |^] (dashed line) and numerical results (circles) 
at a. = 0. Right: relative magnon weight at k = n vs. 8 at 
a = 0. 

In fig. |l|, results are shown for the alternation depen- 
dence of the gap A{5) and the weight of the magnon 
Wpeak(^) at k — 7r relative to the total weight, i.e. peak 
plus continuum. The results at a = are compared to 
the perturbation results of Harris (2^] and numerical data 
of Soos et al. [^l| . The agreement in the region of inter- 
est (S sa 0.03 — 0.06 for CuGeOa) is reassuring to proceed 
along this line. Only for very low S < 0.01 our results 



become incompatible with the the asymptotics [|18 22 
A oc 6 2/3 /y/\ ln(<5)|. The right part of fig. [l] shows the 
spectral weight, which at S = is completely in the di- 
verging continuum Jl5| , p3| . It is shifted very quickly into 
the magnon peak on dimerization. 
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FIG. 2. Dashed curve: magnon dispersion; dotted curve: 
relative magnon weight (right scale); continuum between solid 
lines; dashed-dotted lines: intermediate singularity. Numer- 
ical values: 5 e g = 0.136, J e g = 15.3 meV corresponding to 
5 = 0.0506, J = 11.2 meV. 

In fig. ||, the RPA spectrum is presented. The long- 
dashed line is the magnon dispersion. The magnon 
weight relative to its value at k = n is depicted by the 
dotted line. Note its non-monotonic behavior. The con- 
tinuum is found between the two solid lines. For k < tt/2 
noticeable continuum weight is found only between the 
lower solid line and the dashed-dotted line. The dashed- 
dotted lines show the positions of a singularity between 
the band edges. It results from the local maximum at 
q = tt/2 of the particle-hole energy e p h(g) = e(k/2 + 

q) + e(fc/2 - q) where e{q) = J cS ^J S 2 S sin 2 (q) + cos 2 (g). 
This maximum is not the global one for all k satisfying 

COs(fc) < (l-<5eff)/(l + <5 eff ). 
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FIG. 3. Continua for different wave vectors, the magnon 
5-peaks are not shown. From left to right: upper part 
k/n = 0.01, 0.1, 0.35, 0.5; lower part k/n = 1.0, 0.9, 0.65, 0.5. 

The spectral densities are shown in fig. ||. All singu- 
larities (band edges and intermediate) are of square root 
type. The continua at low k are extremely weak, but ex- 
tend rigorously speaking up to the same high energies as 
do their counterparts at tt — k. For k = (0.5 ± 0.15)7r 
the intermediate singularity (cf. fig. [^) is clearly seen 
at w w 17meV. While for k > tt/2 there is still con- 
siderable weight above the intermediate singularity the 
spectral density for k < tt/2 practically vanishes. The 
extreme asymmetry of the continua at k = tt/2 and 
k = (0.5 ± 0.15)7r can be viewed as the precursor of 
another sharp resonance which separates from the lower 
continuum edge in RPA for slightly larger alternation val- 
ues (S > 0.082). This is a bound state of two magnons. 



The physical origin of this second peak can be under- 
stood best in the limit S — > 1 . Defining A = (1 — §)/ (1 + 5) 
and setting J(l + S) — la perturbation expansion 
around complete dimerization can be performed. To 
assess the influence of the NNN coupling a — a/X is 
kept constant so that A = still corresponds to iso- 
lated dimers. These have k- independent triplet excita- 
tions at uj = 1. The corresponding resolvent R(k,u>) — 
(S z (k)(uj - L)- x S z {k)) is (1 - cos(fc))/(4(w - 1)). Here 
L is the Liouville operator p4|l , which stands for com- 
mutation with H. For A > this mode acquires a 
finite dispersion since the triplets can hop from dimcr 
to dimer with amplitude —1/4. In leading order the 
magnon peak becomes R pca k(k,uj) = po/(uj — u>(k)) 
with 4p = 1 — cos(fc) — A(cos(fc) + 2<5cos(2/c))/4 and 
u(k) = 1 - A(l - 2a) cos(2fcV2. For 5 = 0, Harris cal- 
culated Lo(k) in third order pOf in A. Note that a > 
reduces the triplet hopping. 

A finite A couples the one-triplet subspace also to the 
5=1 two-triplet subspace. This leads in order A 2 to 
the continuum. By a Mori-Zwanzig calculation [53 the 
corresponding spectral density is found in leading order 
to be 

A (l-cos(fc))(l-2acos(fc)) 2 
R cont {k,u) = j 7 === (3) 

where a x = (1 + 2cv)/2, t = (1 — 2a)|cos(fc)| and y is 
the rescaled energy u> = 2 + \y. The spectral density is 
finite for y 2 < t 2 where the square root i n (^) be comes 
iyjt 2 — y 2 . Thus the continuum is cx Xy^t 2 — y 2 /(t 2 + 
a 2 + 2ay). This asymmetric form is qualitatively sim- 
ilar to the results in fig. ||. The edge singularities are 
square roots. Since no phase transition is expected for 
A G (0, 1) this supports the RPA result of square root 
behavior for all finite alternation values in contrast to a 
previous prediction of a divergence fllf . 

The positive value of a\ is the signature of attraction 
between two triplets. A S z = triplet couples to the 
linear combination of two antiparallel S z — 1 and S z = 
— 1 triplets. On adjacent dimers they gain Aa/2 in energy. 
For a\ > t a bound state emerges which is favored by 
a > 0. Its energy is found at the zero of the denominator 
of (^). Since t vanishes at k — tt/2 the bound state is 
most prominent at this value. For a = the bound state 
exists for k € [tt/3, 2-7r/3]; for a > 1/6 it is present at all 
k. 

Similarly one can investigate the two-triplet dynamics 
in the S = 0, 2 subspaces. The secular equation found 
equals the denominator of ([}]) with interactions ai = — ai 
for the quintuplet and clq = 2a,2 for the singlet. This im- 
plies an anti-bound quintuplet above the continuum and 
a more strongly bound singlet below the triplet state for 
all k even at a = 0. The singlet and the quintuplet where 
also found at a = by Harris j2Cj ], The attraction for 
the triplet was missed presumably due to a sign error. It 
is natural to identify this singlet with the -\/3A/2 singlet 
found in the continuum description above. At a — and 



3 



k = 0, 7r the singlet lies at the lower continuum band edge 
which is confirmed numerically |25| for a wide range of 6. 




FIG. 4. Dashed curve: magnon weight relative to total 
weight, right scale; dotted curve: total weight normalized by 
its value at k = tt, left scale. 

Reassured by the coherence of the results of the dif- 
ferent approaches we present in fig. |I] the k dependence 
found in RPA of Wtot(fc), proportional to the equal-time 
structure factor S(k), and of Wp Ca k(fc). The agreement 
with the experimental data [^) is very good. It was pre- 
viously stressed that a « 1/4 is necessary to achieve 
agreement whereas the data in fig. || is computed at 
a = 0. This suggests that S(q) depends mainly on the 
gap A and not so much on a if A is kept constant by 
tuning (5. But the interpretation requires care since the 
approximation tends to raise Wtot(k) around tt/2. Re- 
call in this context the experimental fit result S = 0.042 
which is between our value (5 = 0.0506) and the one used 
so far |8|-|l0|. The fact that for k/n between 0.6 and 0.9 
the spectral weight is found mainly in the continuum is 
important for the experimental determination of S(k) 

Summarizing, we give general arguments for the con- 
tinuum onset at twice the magnon gap in excellent agree- 
ment with experiment Q| . In RPA and for S — ► 1 we find 
a square root singularity at the continuum edge which 
is also in agreement with experiment Detailed re- 
sults for the peak weights and the form of the continua 
are presented. We find evidence in the whole range of 
S > for a singlet excitation below the continuum for 
all k. In addition, we predict triplet (quintuplet) excita- 
tions below (above) the continuum around k = tt/2 for 
not too small S. Their (anti)binding tendency is strongly 
enhanced by competing NNN coupling. An experimen- 
tal identification of either of these excitations would be 
particularly interesting. 

We note that in ref. the continuum is attributed to 
solitonic excitations of the dimerization pattern. How- 
ever, one then has a important deformation of the crystal 
lattice and would expect therefore a considerable increase 
of the effective mass of these objects. The fact that the 



experimental energy scale of the continuum is the bare J 
favors a purely magnetic origin for the continuum. Also, 
in the three-dimensionally ordered structure below Tgp 
solitons on individual chains accompanied by lattice de- 
formations are confined by the interaction with adjacent 
chains and therefore cannot be considered as free excita- 
tions. 
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